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ABSTRACT 
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The increasing interest in long term (inter-annual) 
weather changes and their relation to processes in the ocean 
is beginning to illuminate the need for and the lack of long 
term records of physically significant variables which occur 
over the vast oceanic regions of the northern hemisphere. 

An attempt is made here to evaluate the accuracy of a hind 
cast time series of surface wind vector components and speeds 
for the period of January 1960 through December 1969. 

The tano hic model used to calculate these 
records is described as well as the twelve-hourly surface 
pressure data which were used as input. The central moments 
of the probability distributions of the computed records are 
compared to those of corresponding time series observed on 
Ocean Station Vessels. Linear correlation coefficients 
between observed and computed records were found to average 
0.8] for the components of the wind vector and 0.65 for the 
wind speed. Regression relations between computed and 
observed records also are presented. Spectral analysis of 
low pass filtered records showed coherence between observed 
and computed records increased with decreasing frequency. 
The accuracy of the computed records in low latitudes also is 


estimated. 
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I INTRODUCTION 


A. GENERAL 

One of the most important variables encountered in the 
study of physical oceanography and air-sea interaction is 
the wind field at the sea surface. To a large degree, the 
surface wind velocity and its variability, both in time and 
Space, determine the exchange of heat and momentum across 
the air-sea interface. Through such exchange the wind affects 
not only oceanic surface phenomena, such as waves and surface 
currents, but also features within the sea such as the 
vertical thermal structure and deep currents. 

The smallest scale of concern in this study is the 
Synoptic scale. Such a scale covers a period of time from 
several hours to several days and a distance from several 
tens of kilometers to several hundred kilometers. Past 
attempts at estimating the surface wind field on this scale 
have been of two types. One could either average wind obser- 
vations collected over a given area for a specified period 
of time (Seckel, 1970) or one could estimate the surface winds 
from a similar collection of observations of other variables, 
such as surface atmospheric pressure (Roden, 1974). 

The accuracy of synoptic scale surface wind fields 
derived directly from wind observations suffers from two 
major drawbacks: a paucity of observations, and the high 


variability of the small scale wind field. The small number 
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of observations usually available in a typical mid-ocean 
region over a three-or-four-day period are very often grouped 
together in a small section of the region, often along major 
shipping lanes or along the track of a single ship. In 
addition, major portions of the open ocean lack any reports 
of wind conditions for weeks at a T As a result, the 

few reports available must necessarily be used to describe 
large areas when estimating the synoptic scale wind field. 
Verploegh (1967) found, however, that simultaneous wind speed 
observations from ships less than sixty miles apart in 
"synoptically homogeneous areas" differed significantly. His 
analysis indicated the correlation coefficient between such 
observations was less than 0.7. 

Because of such considerations, the second, indirect 
approach is generally used to estimate the synoptic scale 
surface wind field. This approach depends on the fact that 
the distribution of surface atmospheric pressure can be more 
accurately described than can the surface wind field from a 
limited number of observations, due to the lower variability 
in time and space of the surface atmospheric pressure field. 
An accurate estimate of the synoptic-scale surface wind field 
may be made by applying the geostrophic approximation to the 


surface atmospheric pressure Field. 
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B. OBJECTIVES 

Although the relation of the geostrophic to observed 
wind is a popular topic, this study is not intended to be a 
detailed examination of that relation. Rather, it is limited 
to an evaluation of one particular surface wind model used, 
in this study, to hindcast time series of surface winds over 
oceanic areas. The model is quasi-geostrophic, incorporating 
the effects of friction. It is one part of operational 
oceanographic models used at Fleet Numerical Weather Central 
and at Fleet Weather Central, Rota (Spain).! 

As a result of this study, the ten-year average and 
linear trend of observed wind conditions at several points 
in the North Atlantic and North Pacific Oceans were determined 
and compared to those calculated from the model. The primary 
ect veror the study was however, to determine the accuracy 
of the computed time series and to determine if calibration 
factors could be applied directly to them in order to increase 


their accuracy. 


‘The model was developed by the author while at Fleet 
Numerical Weather Central; most of the work of this study was 
carried out while the author was at the Environmental Prediction 
Research Facility. 





II. APPROACH 


A. THEORETICAL BACKGROUND 


E 


The Geostrophic Mind Equation 


Geostrophic wind has been considered to be an 


acceptable estimate of the true surface wind in many studies 


of oceanography and air-sea interaction (Roden op. cit., 


Namias, 1963). Geostrophy assumes a balance between the 


pressure gradient force and the Coriolis force and ignores 


acceleration, friction, and vertical motion. This balance is 


given in equation (1: 


h V 
wnere g 
p 
EH 

X 


H 


k 


T 
we k X -7 VP (1) 


geostrophic wind vector 


surface pressure 


CN MEUS 
9X ay 


vector product operator 


j (horizontal gradient operator) 


20 sin @ (the Coriolis parameter) 
angular rotation of the earth 
latitude 

air density 


unit vertical vector 


Although this approximation is widely used, its accuracy is 


limited by the required assumptions indicated above. 





The true surface wind vector is the resultant of 
the geostrophic wind vector and an ageostrophic wind vector. 
The direction and magnitude of the ageostrophic wind vector 
determines the accuracy of the geostrophic approximation. 
Unfortunately the ageostrophic vector has been shown to vary 
significantly over the Northern Hemisphere (Roll, 1965; 
Carstensen, 1967; Brümmer, et al., 1974). Such variability 
makes the geostrophic wind hard to relate to the surface wind 
in studies on synoptic or larger scales. 

The accuracy of the geostrophic approximation is 
also limited, in practical application, by the temporal 
variability of the surface pressure field. In practice, the 
geostrophic approximation is often applied to surface pressure 
distributions averaged over a month or more (Namias, op cit.). 
Such a procedure results in calculating the vector resultant 
of the geostrophic wind for the averaging period. Obviously, 
non-conservative phenomena which depend on the surface wind 
cannot be treated using such a resultant wind field. In the 
study of such phenomena, the resultant vector is assumed to 
represent the mean or steady-state wind prevailing over the 
averaging interval. 

Whether or not the resultant wind is a good esti- 
mate of the mean wind depends on the variance of the wind 
vector over that interval. The mean wind speed, Wyo and the 


resultant wind speed, Wp are given below in terms of components: 
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^ 
PO e£) 
Wy = 2 (2) 
t=] T 
"RU PI 
WEINE De = (3) 
R ee tid 
where gd? V. are orthogonal components of the wind vector at 


time t for 1 « t « T. (Averaging twelve hourly data over 


thirty days, we have 1 « t « 60). The difference between 


We and We can be expressed in terms of the variance of each 


of the components, of and st , and of the wind speed, of; 


2 2 2 2 2 
Wy - Wp E + EE (4) 
, tu ft u,\? 
2 _ t EE ) 
where cM M rs Lo e (5) 
t=] t=] 


and similarly for of, and st. As will be shown in Section 


III-A, the difference expressed in equation (4) results in 
a drastic under-estimate of wind speed when the resultant 
wind is used in place of the mean wind speed in problems 
involving the square of the wind speed. 

2 The Surface Mind-Geostrophic Mind Relation 

Although this is not a study of the general relation 

of the geostrophic wind to the surface wind, a few words 
regarding this relation are appropriate for background. As 


indicated, the surface wind vector can be resolved into 


(ES 





geostrophic and ageostrophic components. The ageostrophic 
component can be assumed to consist of a S Oh component 
and an acceleration component (Brümmer, et al., op cit.). 

The frictional component is generally accepted as 
the dominant ageostrophic component (Verploegh, op cit., 
Brümmer, et al., op cit.). Haltiner and Martin (1957) suggest 
that the frictional component is dependent on the stability 
of the surface layer of the atmosphere. Hasse (1974) has 
demonstrated that the dependence of the surface wind speed on 
the surface pressure gradient is an order of magnitude greater 
than its dependence on stability. 

The effects of stability have not been explicitly 
incorporated into the model; to some degree however, they 
have been implicitly incorporated through the latitude 
dependence of terms in the various equations of the model. 

The major effects of stability are assumed to occur in the 
region between 45°N and 25°N. This implicit incorporation of 


Stability is shown in the ratio of the surface wind speed to 


the geostrophic wind speed in Figure 1. 


B. METHOD OF COMPUTING THE SURFACE WIND 

The calculations were performed on the Fleet Numerical 
Weather Central (FNWC) 63x63 polar stereographic grid of the 
northern hemisphere (FNWC, 1974). On this square grid the 


equator is an inscribed circle, i.e., tangent to the grid 





boundary. The grid mesh length at 60°N is 38] km and increases 
with increasing latitude. The data, although calculated on a 
polar stereographic grid, will be discussed relative to a 
Mercator grid.“ 

The components of the horizontal pressure gradient at a 
given grid point (i,j), were approximated to the second order 


by equation (6): 


Zei? A! = BECA KıtPli-2,j) - K,LP(i-1,5) 


(6) 
SPUR ees 
where Ky = 1/12 
K, =8 
L - 381 km (mesh length) 


LESIN?60°) 


MF(i,j) = (STE - map factor of the grid at 
point (i,j) 
D = latitude of grid point (i,j), 


and similarly for 271,3). Using (1) with a constant air 


3 


density of 1.22 x 10 gm cm ^, the geostrophic wind vector 


“Both Mercator and polar stereographic projections are 
conformal projections of the earth and as a result, angles 
are preserved when transformed from one to the other (Taylor, 
1955). 
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3 The magnitude of this vector (the 


at each point was found. 
geostrophic wind speed) was then reduced as a function of the 


latitude of the grid point according to (7): 


^ 


VG 3) 5 A EB(9)] Vo CH L3) (7) 
where NEE = surface wind speed at grid point (i,j) 

Vo O3) 
A = 0.93 


geostrophic wind speed at grid point (i,j) 
4 


(0.65+0.2(9/25°); 9 < 25° 
B(g) IT ; 259 « 9 « 45? 
0.75+0. 1208), 450 <p 


The resulting ratio of (VIV) > which is strongly 
dependent on latitude, was suggested by data from Rol] 
(op cit.) and Carstensen (op cit.). The value of NEEN 


used in the model is plotted in Figure 1 along with values 


The sine function, SILA, used in the calculation of the 
Coriolis parameter for the geostrophic wind calculation and 
in the calculation of the deflection angle (equations ] and 
8) was modified below 35?N. This was done in an attempt to 
reduce the impact of this parameter in the calculation of 
equatorial winds. The modification was suggested by the 
subroutine SNLTSRX in the Fleet Numerical Weather Central 
Prodr amali brary. (FHWC. op cit.) ihe value of SILA south. of 
35°N is given by the tangent to the sine curve at 35?N, i.e., 
SILA(g) = 0.0144 (SIN-1(0]) NIU AEDEM Pe 35". 


ee factor A resulted from erroneously using 408 km 
for the grid mesh length in the model (D. Dale, FNMC, 
personal communication). 
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found by Roll and by Carstensen. Although Carstensen derived 
his relation from a much smaller data set than that in Roll's 
study, the major features of the latitude dependence are 
similar. The value of SI used in the present model 
attempts to reproduce the first order dependence on latitude: 
i.e., the maximum in the mid-latitudes, sharply decreasing 
with decreasing latitude toward the equator and a more modest 
decrease with increasing latitude above 45°N. 

After being reduced in speed, the surface wind vector 
is then rotated to the left of the geostrophic wind vector 
(toward low pressure). Together, deflection and reduction 
of the geostrophic wind vector simulate the effects of friction 
(Haltiner and Martin, op cit.). The deflection angle, o, was 


calculated from (8): 


2 


a = Ka (K, - KN. / (1+SIN Q) (8) 
where K4 = 1.475 
Ko = 22.5 
E -2 
K3 AS | 0 
V, 7 surface wind speed, from (7) 


Figure 2 is a plot of a against latitude for several values 
of surface wind speed. 

The dependence of a on the square of the surface wind 
speed strongly decreases the deflection angle associated with 
high surface wind speeds. This approach was based on the 


assumption that high winds (associated with well-defined 
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cyclones) generally show less cross-isobar flow than do the 
low winds associated with nearly flat surface pressure 
distributions. The assumed reduction in the deflection angle 
as wind speed increases could be caused by a slight reduction 
in the frictional drag experienced at higher values of surface 
wind speed. In this model, however, this reduction becomes 
significant only at relatively high wind speed. For example, 


] 


at 20 m sec the deflection angle is 75% of its value at 


] 


1 m sec^!. At 30 m sec! the deflection angle has been 


reduced to 30% of the value at Im sec”. 

As with the value of SAPE the latitude dependence 
of @ also was suggested by Roll and Carstensen. Figure 3 is 
a plot of the value of a against latitude as found by Roll, 
by Carstensen, and as used in this model. The values of « 
from the model are those for surface wind speeds of 5 and 
15 m sec”! Less than 15% of the reported wind speeds in 
this study were outside this range. As in the case of the 
reduction factor, the expression for a attempted to reproduce 


the major features of the latitude dependence shown by both 


Rolà's and Carstensen's data. 


Cz DATA GENERATION 
E Input Data for the Model 
The model was used to calculate a ten-year time 
series of the surface winds at twelve hour intervals at grid 
points near international Ocean Stations. The input for the 


mode] consisted of surface pressure fields on the FNWC grid. 
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A11 the available surface pressure analyses for synoptic 
periods of 0000 GMT and 1200 GMT were used for the ten-year 
period from 0000 GMT 1 January 1960 to 1200 GMT 31 December 
1969. However, not all the analyses for this period were 
available and this required some conditioning of the output 
data which will be discussed in Section II-D. 

The number of surface pressure analyses used for 
each month in the ten-year period are shown in Table I. In 
general, the analyses were fairly evenly distributed over 
the twelve months. The month of April had the poorest 
representation with 488 (812) of the 600 possible analyses 
used, while November was the best represented with 534 (897) 
of the 600 possible analyses. Most of the missing data were 
concentrated in the period from April 1960 through June 1962. 
Within this period only one analysis per day was available. 
As a result, only 365 (50%) of the 730 possible analyses were 
used for 1961. The best-represented year was 1967 for which 
728 (99.72) of the 730 possible analyses were available. The 
missing analyses totaled 1098 (15%) of the 7306 possible in 
the ten-year interval. 

According to Bakun (1973), who used the same set 
of surface pressure analyses to calculate coastal upwelling 
indices, the analyses between January 1960 and July 1962 
WerpEomnetuae National Climatic Center at Asheville, North 
Carolina. The analyses from July 1962 through December 1969 


were produced by objective computer analyses at FNWC. 
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2 The Observational Data 

The observed wind data were taken at six interna- 
tional Ocean Stations (0S). These were the only continuous 
time series of observed wind speed and direction available 
for the ten years of this study. As these data were assumed 
to be of high quality, no attempt was mae to determine their 
accuracy other than gross error checking as described in 
Section II-D. The observations were thus used as the standard 
against which the calculated wind records were verified. 

The observed wind speed and direction were extracted 
from the Surface Marine Observation Tape Family 11 from the 
National Weather Record Center for the following Ocean 
Stations: BRAVO (OSB), DELTA (OSD), JULIET (0SJ), PAPA (OSP), 
VICTOR (OSV), and NOVEMBER (OSN).° Table II indicates the 
position of each OS, the position of the computational grid 
points, and the distance between each OS and its correspond- 
ing grid point. 

SH Data Used in Analysis ` 

Both the observed and computed data were transformed 
to zonal (East-Mest) and meridional (North-South) components. 
The resulting set of three time series for each OS and each 


grid point (a total of 39) consisted of the values of zonal 


?A11 the Ocean Station records with the exception of 
NOVEMBER (OSN) were compared to computational time series 
irom a single grid point near their location, Time series 
were calculated at two grid points bracketing OSN in order 
to estimate the variability of the computed wind over a 
single grid mesh length. These computation grid points will 
Be relerred to, henceforth, as NIC and N2C. 
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and meridional components and the wind speed at 0000 GMT and 
1200 GMT for the ten-year period from | January 1960 to 31 
December 1969. It should be noted here that a negative zonal 
value indicates a wind from the West and a negative meridional 


value indicates a wind from the South. 


D. DATA CONDITIONING 
l. Gaps and Gross Errors 

A check for gross errors was performed following 
extraction and transformation of the observed and computed 
wind vectors. A wind speed of 50 m sec! was selected as the 
maximum allowable value; this established a gross error 
criterion in order to eliminate only those errors caused by 
hardware failure such as tape parity. No wind speed values 


| were found in either the computed or 


in excess of 50 m sec 
observed data. 

After gross error checking was completed, ensemble 
analysis was performed on the observed and computed time 
series. Following that, the gaps in the records were removed. 
The computed records differed markedly from the observational 
data in both the number and length of gaps. For the eighteen 
observational records, the number of single point gaps per 
record averaged about 200, while about 300 points per record 
were missing in gaps of two or more consecutive points. Only 


about 150 points per record were missing from the computed 


time series in the form of multiple point gaps. However, for 
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these computed records there were about 950 single point 
gaps, due primarily to missing pressure analyses within the 
period of April 1960 through June 1962. 

Because of the high proportion of single missing 
data points in the calculated record, a single missing value 
was replaced by linear interpolation between the preceding 
and following value. For gaps consisting of two or more data 
points, the missing data were replaced by the time series 
average. This procedure was suggested as a means of filling 
gaps while not introducing spurious data into the records 
which would be subjected to subsequent spectral analysis.° 

23 Filtering 

Filtering was done on the computed and observed 
records in order to allow comparison of records with near]y 
equal data content. Because of the substantial missing data 
in the high frequency part of the computed record (due to 
single point gaps) and because of the stated interest in 
large scale phenomena, a low-pass filter was applied to each 


of the computed and observed time series after the gaps were 


SE. Thornton, Dept. of Oceanography, Naval Postgraduate 
School, personal communication. 
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removed. A five-point binomial filter indicated in equation 


(9) was used.’ 
2 
Meo) - F a. X(ttiAt) (9) 
12-2 
where Y(t) z filtered value at time t 
a; = binomial coefficient (0.38, 0.25, 0.06, 
A A 81,42) 
X(t+iAat) = unfiltered value of the time series at 
time t+iat 
At = sample interval (0.5 days) 


Figure 4 shows the effect of applying this filter to a typical 
segment of one of the computed time series. The amplitude of 
the residual series (unfiltered minus filtered) is approxi- 
mately the same amplitude as the filtered series, indicating 
that approximately 50% of the variance of the unfiltered 

series was removed by filtering. This estimate of 50% was 
confirmed by a comparison of the filtered and unfiltered record 
variance (not shown). Before any E or gap removal was 
done, however, the observed and computed records were analyzed 


as ensembles in order to determine the relation between them. 


such a filter has a frequency response, R(f), at 
frequency f 4 
Pit) paces (iat) (10) 
(Panofsky and Brier, 1958). 
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IU RESULTS 


A. ENSEMBLE ANALYSIS 

In this section the ensembles of the observed and 
computed data are compared. The term ensemble is taken to 
mean the set of values which constitute a given time series 
of the variable in question. Due to the length of the records 
(ten years), significant questions about short-term accuracy 
of the model cannot be treated by ensemble analysis. The 
problem of short-term accuracy is investigated through 
spectral analysis as reported in Section III-B. 

ie Central Moments 

Doth the observed and computed data were considered 

to have come from specific probability densities. While these 
probability densities were not explicitly defined in this 
study, they were described by determination of their central 
moments. The general form of the central moment of a given 


probability density is given by equation (11). 


GZ Ca" GER 
RENE CH 
where up 7 pin central moment 
eo ¡En value of the variable X 
N = total number of values of X in the ensemble 
uc x,/il 


i=] 


= mean value of the ensemble 
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In practice, the third and fourth moments are 
normalized by the third and fourth powers of a (a = Uy) > 
the standard deviation of NE ensemble (Tennekes and Lumley, 
1972). In this study the first moment is taken to be the 
mean rather than zero as in equation (11). 

The criterion used to estimate the accuracy of the 
computed data - the "figure of merit" - was the ratio of a 
central moment of the computed record to that of the observed 
record. This simple criterion was intended to show only two 
things: the percentage of a given observed statistic present 
in the corresponding computed record; and those statistics 
which differed in sign between computed and cbserved records 
(resulting in a negative figure of merit). 

As the computed data were compared to observations 


8 os, the average figure of merit for the model was 


at six 
taken to be the average of the figures of merit for each OS. 
An average was found also for the three OS in the Pacific and 
for the three in the Atlantic. 

a. Wind Speed 


The values of the central moments and their 


figures of merit for the observed and computed wind speed are 


The two computed records, NC1 and NC2, were compared 
to the observational record at OSil. In order not to over- 
emphasize these data when considering averages, the average 
of the figures of merit of NC2 and NC] was used as the figure 
of merit at OSN when calculating the overall average figure 
of merit or the average of the figures of merit in the Pacific. 





tabulated in Table III and are plotted against latitude in 
Figure 5. The average computed wind speed (first moment) 
varied between 95% and 91% of the average wind speed observed 
at the six OS. Over the latitude range of the 0S (30°N to 


56.5°N), the average observed wind speed decreased from a 


] ] 


maximum of 10.6 m sec . in the Mesterlies (OSP) to 6.6 m sec. 
in the North East Trades (OSN) for the ten years between 1960 
and 1969. The average figure of merit for the mean computed 
wind speed is 0.92. This figure is 0.9] in the Pacific and 
0.93 in the Atlantic. 

It is worth noting again here that we are 
dealing with an average wind speed and not a resultant wind 
speed. The resultant wind speed can be determined by 
combining tne average zonal and meridional components by 
vector addition. The average wind speed was determined from 
the wind speed at each interval in the ten year time series. 

The variance (second moment) of the wind speed 


e er in the latitudes of 


the Westerlies to approximately 10 n* Bc in the latitudes 


decreased from approximately 28 m 


of the Trades - a latitude dependence similar to that of the 
average wind speed. The average figure of merit for the 
Damıance is 1.072. Only at 1s the figure of merit (0.89) 
less than unity. The maximum (1.24) is shared by OSJ and OSN 
(NC2). Mhen comparing the variances of the computed and 
observed records by means of the F test it becomes apparent 


that any figure of merit for the variance which is not equal 
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to unity is statistically significant at the 12 level due to 
the very large (56000) number of degrees of freedom for each 
ensemble (Freund, 1962). That is, the probability that the 
computed and observed ensembles come from the same probability 
density is less than 0.01. Relating these statistically 
Significant differences in the Ee eerie computed records to 
to physical causes is, for the most part, beyond the scope of 
this study. However, the location of the grid point NC2 

(at 32N) at the northern boundary of the North East Trades 
appears associated with the high figure of merit at this 
location. Seckel (op cit.) suggested the wind speed is much 
more variable at 32N (NC2) than it is in the Trades (at OSN), 
and thus a figure of merit for wind speed variance greater 

than 1.0 is to be expected between these records. This suggests 
that the model can resolve significant changes in the varia- 
bility of the wind field which occur on a scale of one 
grid-mesh length. 

The standard deviation (square root of the 
variance) can be used to estimate an upper bound on the 
average fluctuation of the wind speed from the ten-year mean. 
According to Chebyshev's Inequality, the probability that the 
wind speed will differ from the average by more than 1.414 
standard deviations is less than 0.50 (Freund, op cit.). The 
value of 1.414 o in the Westerlies was about 7 m 3989 This 


value is considerably greater than the diurnal range of the 
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wind speed in the open ocean ODORE Diy RO OP cite): 
Thus, most of the variance in the wind speed observed at OS 
over the ten-year period was due to variations of a period 
longer than twelve hours. The 505 reduction in the variance 
due to filtering as indicated in Section II-D-2 suggests that 
50% of the variance is due to wind E. fluctuations of a 
period greater than two days. 

The third moment (skewness) of the computed 
wind speed is higher than that of the reported wind speeds 
at all the OS. The average figure of merit for all the 
locations is 1.30. This high value represents a tendency of 
the computed wind speed to exceed reported wind speed 
Significantly in cases of high wind. Of the locations 
studied, this tendency is most pronounced in the Pacific 
Trades. The fact that the average figure of merit was greater 
than one for the skewness and less than one for the average 
of the wind speed suggests that the overestimated high wind 
speeds in the computed record were more than compensated by 
more numerous underestimations of low wind speeds. The 
omission of the effects of curvature, i.e., the acceleration 
of the wind due to the centrifugal force, in the geostrophic 


equation was suggested as the cause of this effect." 


?R. Elsberry, Dept. of Meteorology, Haval Postgraduate 
School, personal communication. 


de 





The exaggeration of both extremely low and 
high wind speeds is also suggested by the fourth moment 
(kurtosis) which indicates the amount of extreme data in the 
ensemble, i.e., the number of data points in the ensemble 
which are far from the mean. The figure of merit for 
kurtosis is highest for the Pacific in the transition from 
Westerlies to Trades (OSV and OSN-NC2). As can be seen in 
Figure 5, the value of the kurtosis of the wind speed 
increases dramatically in the Trades. However, this increase 
is an artifact of the normalization by the fourth power of 
the standard deviation, which is much smaller in the Trades 
than in the Westerlies, and does not indicate that very strong 
winds occur more frequently in the Trades than in the Westerlies. 

D. Zonal Component 

The first four moments for the zonal component 
of the wind vector are shown plotted against latitude in 
Figure 6. Their values and figures of merit are tabulated 
in Table IV. The latitude dependence of the first moment 
(average or resultant) of the observed and computed zonal 
components is similar. The difference between the first 
moments of the observed and computed zonal components is 
greatest at OSP (0.85 m sem dy While the average figure of 
Herat of the first moment is 1.0, the model is more accurate 
in Age Atlantic than in the Pacific Ocean. In the Pacific the 
figure of merit ranges from 0.43 to 1.45 and averages 1.01. 
Inm Mike Atlantic the range was smaller, from 0.89 to 1.05, 


and the figure of merit averages 0.98. 


36 








This disparity between Atlantic and Pacific 
is evident also when comparing the ae (or second moment) 
of the observed zonal component to that of the computed zonal 
component. The figure of merit ranges between 0.83 and 0.95 
in the Atlantic Ocean and between 0.78 and 1.02 in the 
Pacific. Overall, 1t averages 0.88, indicating that 
the model reproduced about 88% of the variance of the zonal | 
wind. The relatively poor accuracy of the model at OSP as 
shown in Figure 6 is unexplained; however, the poor agreement 
between OSN and grid point NC2 is probably due to the transi- 
tion to the Trades as previously suggested. 

The skewness showed the poorest agreement of 
all the moments calculated for the computed and observed 
zonal components. The model consistently underestimated this 
parameter; the figure of merit averaged only 0.63. The signs 
of the computed and observed skewnesses differed at QSJ; 
however, the values at this point were so near zero that the 
difference is insignificant. 

The kurtosis (fourth moment) of the computed 
zonal component was consistently higher than that of the 
observed zonal component. The average figure of merit was 
ies andewas nearly equal in both the Atlantic and Pacific 
Oceans. 

e Meridional Component 

The first four moments of the probability 
density of the computed and observed een components 
are gabulated in lable V and are plotted against latitude in 


Figure» /. 
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The figure of merit for the mean of the 
meridional component, averaged for all locations, was 1.22 
and was nearly equal in the Pacific and the Atlantic Oceans. 
The individual figures of merit of the means suggest that the 
computed meridional component has a bias of 0.3 to 0.9 m sec’!. 
This bias is a fictitious (computed) component from the south, 
which reduced the average computed meridionaj component at 
those stations where the resultant direction was from the 
north (0SB and 0SN); and it increased the meridional component 
at those stations where the resultant direction was from the 
south. Such a consistent bias suggests that a fictitious 
east-west gradient (higher pressure to the east) in the sur- 
face pressure analyses may occur over both the Pacific and 
Atlantic Oceans. 

The variance of the observed and computed 
meridional components showed only slightly better agreement 
than was the case for the zonal component. Averaged over the 
locations studied, the model reproduced 92% of the variance 
of the observed meridional component. Again, OSN (NC2) 
showed the highest figure of merit due to the effects of the 
transition to the Trades. 

The skewness of both the observed and computed 
values was nearly zero at all locations north of the Irades. 
This indicates that, relative to the average meridional 
component, the strong positive (northerly) and negative 


(southerly) fluctuations were equally distributed. The 


35 





negative value of the skewness at OSN indicates that the 
southerly (relative to the mean) extremes predominated. 
Because the skewness was so near zero for all stations north 
of the Trades, figures of merit at these locations have little 
meaning. In the Trades, the model overestimated the skewness 
by up to 98%. This very high figure na is ın marked 
contrast to that for the skewness of the zonal component 
(Table IV). This discrepancy suggests that in the Trades, 
extreme values in the calculated record of meridional com- 
ponents are not always matched in the observed record. 

The value of the kurtosis (fourth moment) of 
the meridional component was slightly overestimated by the 
model. The figure of merit averaged 1.19, and the averages 
for the Atlantic and Pacific Oceans were identical. 

The values of the skewness and kurtosis of the 
meridional component north of the Trades, approximately O and 
3.0 respectively, suggest that the values of the meridional 
component exhibit the well known Gaussian distribution in 
enS area. 

Of the four central moments discussed above, 
perhaps the most important for oceanographic purposes is the 
variance. This parameter is a measure of the kinetic energy 
Othe wind, and is directly related to the frictional 
coupling of the ocean and atmosphere. 

On the average, for the locations studied, the 
model reproduced about 90% of the variance in each of the 


components of the wind vector and overestimated the variance 
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of the wind speed by 7%. In none of the cases was the error 
greater than 25%. The consistent overestimation of the 
variance at NC2 as compared to NCI indicates that features 
such as the transition from the Trades to Westerlies can be 
resolved within one mesh length of the computation grid. 
The difference, shown 2 GUN o 4, between 
the square of the mean wind speed (derived from individual 
pressure analyses as in this study) and the square of the 
resultant wind speed (derived from monthly mean pressure 
analyses) is often neglected in studies of large scale air- 
sea interaction (Namias, op cit.). As can be seen from 
Tables III through V, this difference as observed at OSP jis 


92.5 n* ae The resultant wind squared ıs only 


20.7 m^ sec ^. In this case the resultant square is less 


than 20% of the mean squared (112.2 n* dE 
The difference between the ten-year mean and 
ten-year resultant is of course larger than the corresponding 
difference over a month due to the longer period (seasonal 
and year-to-year) variation. However, this example does show 
that the difference between the square of the mean and the 
square of the resultant can not be considered insignificant 
when dealing with the energetics of large scale air-sea 
interaction. 
2: Ten-Year Trend 
The estimation of a time series of wind conditions 


using pressure analyses from different sources raises a 


problem cited by Bakun (op cit.). That is, any long-term 





trend in the wind record may be due to different methods of 
analyzing the pressure distribution in successive time 
intervals. For example, if improved analysis techniques or 
more numerous pressure reports allowed better resolution of 
tight pressure gradients there could be a fictitious increase 
in the wind speed derived from the improved pressure analyses. 
Obviously, such a fictitious increase should not be reflected 
in the long-term trend of observed wind speeds. 

In order to investigate this possibility, the 
ten-year changes in the observed and computed wind speeds and 
components were calculated using equation (12) (Bendat and 


Persol, 1971). 


u(t) s b. + by t (12) 
U 
where u(t) = value of trend line at time E 10 and 
N H. 
2(2N+1) zu - 6 znu 
nel " E 
Dr = (13) 
N(N-1) 
N N 
12 E nu - 6(N+1) Zu 
Ba EHS 
Bue (14) 


h N(N-1) (N+1) 


is sample interval, h, in equation (14) is assumed 
Fo be equal to |, However, the gaps in the record were 
not filled prior to the determination of the ten-year trend. 
Because the vast majority of the gaps consisted of a single 
missing data point, this should have little impact in the 
determination of the ten-year change. 
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uc nen value of the ensemble 

N = total number of values in the ensemble. 
The change over the ten-year period, A, is given by equation 
(15): 


Ni (T5) 


The values of A for the wind speed and each of the components 
are shown plotted against latitude in Figure 8. 

In the Westerlies the reported wind speed has 
decreased over the ten years between 1960 and 1970. In the 
Trades (OSN) it has increased very slightly. The decrease 
appears to be most pronounced in the Atlantic Ocean with a 


l at OSJ. For the other 


] 


maximum decrease of 1.16 m sec 
locations the ten-year change was less than 0.6 m sec. 
Computed wind speeds showed changes similar to those observed 
at all locations with the exception of OSP, where the sliaht 
increase of the computed speed (in opposition to the observed 
decrease) is unexplained. 

With the exception of OSJ and OSN the values of the 
observed zonal components show a general increase over the 
ten years of the study. It should be noted that, due to the 
negative sign of the prevailing zonal comporent in the 
Westerlies and the positive sign in the Trades, this repre- 
sents a decrease in the prevailing zonal components at all 


tnewibcationswe»xcept at 059. 
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At OSJ the observed increase in the prevailing 
zonal component was not accompanied by an increase in the 
prevailing meridional component. In fact, the prevailing 
observed meridional components decreased at all locations, 
with the exception of OSD. 

Although there are noticeable differences between 
the trends of the computed and observed series, i.e., OSV 
(zonal), OSD (meridional) and OSP (speed), the similarity in 
latitude dependence and magnitude between observed and 
computed records suggests that the trends in the computed 
series are largely controlled by processes controlling the 
trend of the observed series, and are not due to artificial 
changes in the pressure data. 

oF: Correlation and Regression Analysis 

This section deals with the linear correlation and 
non-linear regression analyses carried out on the ensembles 
of paired (computed and observed) data.  Non-linear 
(parabolic) regression was done in order to determine the 
non-linearity of the relation between the computed wind data 
and those observed at the OS. 

The linear correlation coefficient, R, is found 


using equation (16): 
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where SE computed value 


yn dm observed value 


N X 
- n 
x = 2 (MW 
n=] N 
E N y 
M E UL Ee 
n= ] 


The value of Re 1S a measure of the dependence of the varia- 
bility of one variable on the variability of the other through 
linear coupling (Panofsky and Brier, op cit.). Table VI 
indicates the values of R and R? for the wind speed and 
components, and the number of paired values upon which the 
correlation is based. In all cases, the correlation between 
computed and observed values of the components exceeds that 
between computed and observed speed. Averaged over the six 
OS, the model explained 65% of the variance in each of the 
components and 42% of the variance in the wind speed. 

The lower correlation for wind speed can be 
attributed to the combined effect of errors in both the zonal 
and meridional components. If the computed wind vector is 


^ 


W., then 
C 


Ez NM S (17) 


where Le and A are the computed zonal and meridional con- 
ponents and i, j are unit vectors in zonal and meridional 


directions. If we also define the wind error vector M ds 
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W = W - W 3 (18) 


where M. is the observed wind vector 


Wo = Zi + a x (19) 


we have, in terms of components, 


D zZ) Ir mm) (20) 


C C 


and the wind speed error (magnitude of the error vector) is 


given by 


» 
+ (M_ - M )*) (209 


C 
Thus, the error (and resulting lack of correlation) in the 
wind speeds should always be greater than or equal to the 
error (and lack of correlation) in the components; i.e., 

the correlation coefficient for the wind speed should always 
be less than that for either of the components: 

On the average, the computed data correlated best 
with observations for the zonal components. In the Atlantic 
Ocean 67% of the observed variance in the zonal component was 
explained by the model. However, the meridional component 
Showed higher correlation than the zonal component for five 


of the seven series. 
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In addition to linear correlation, non-linear 
regression between the computed and observed records were 


also determined. The regression equations are given below: 


a 2 
0 1 


-£ 
i! 


X * AX (22) 


B + B.Y +BY? 


0 1 2 (23) 


>< 
Hu 


where Y = computed value 


observed value 


>< 
I 


The coefficients were found by solving the follow- 
ing set of simultaneous equations (24a - 24c) (Spiegel, 1961) 
for the computed data as a function of the observations, and 
solving the inverse set (replacing Y by X and solving for 
B, B. and B 


: >) for the observations as a function of the 


computed data. 


2 


Pu = AN + Ay EX + A EX (24a) 
d Z 3 
ENY E A eek NX (24b) 
O 1 2 
EXÓY = AQEXÓ € AQrX? Ap (24c) 


Table VII lists the coefficients Ao? Ai, and A, for equation 
(22). Table VIII gives the coefficients for the inverse 
equation (23). The values of N, the number of paired data 


values in each series, are found in Table VI. 
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Comparing A, to A, and B} to B, indicates that an 
essentially linear relation.exists between the computed 
records and observed records. The non-linear coefficients, 
A, and Ba > are two orders of magnitude less than d and B 
for all series except for the meridional component at OSN. 
At this location the non-linear term was only one order of 
magnitude less than the linear term. 

Figures 9 and 10 are plots of the two regression 
equations for wind speed at OSJ and OSN (NC2). These loca- 
tions showed the highest and lowest correlation between 
observed and computed wind speed. The separation of the two 
regression curves is due to the scatter of the computed and 
observed wind speed data being compared, as shown in Figures 
11 and 12 (together with Table IX). The scatter of the data 
is at least partially due to neglecting some of the signifi- 
cant physics in the model, i.e., stability, and the curvature 
of the isobars in the pressure analyses. The area where the 
regression lines converge (the area of least scatter) is the 
area of best fit of the model. For OSJ this region is 


centered around 11] m see; for OSN it centers around 


7] m sec. In all cases this region of best fit is near the 
mean value of the series. 

If the relation of the computed parameter to the 
observed parameter changes with time, then no single time- 
independent equation can express this relation. For the 
purposes of describing the relative accuracy of the model at 


various locations, however, the relation between the 
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observed and computed data is assumed to be shown by a single 
"true" regression line. The "true" regression line is defined 
as that line bisecting the distance between the two regression 
curves X = f(Y) and Y = f(X). This distance is measured per- 
pendicular to the X=Y line. The "true" regression lines were 
found by graphical means for the wind speed at OSJ and OSN 
and are plotted in Figures 9 and 10. 

These figures indicate that a consistent feature 
of the model, overestimation of high observed wind speeds 
and underestimation of low values, is most pronounced in 
the low latitudes (OSN) where observed wind speeds below 
about ]l m sec” are underestimated. This underestimation 
occurs to a lesser degree at OSJ up to wind speeds of about 
14 m sec”). The scatter cannot be reduced by applying 
correction factors (deduced from the regression analysis) to 
the output of the model; however, taking average values of 
the model output over periods longer than 12 hours will 
provide a more accurate time series (at the expense of 
reduced time resolution). This is a result of the elimination 
of the extreme values in the time series by the averaging 
process. As the averaging interval increases, the mean over 
that interval approaches the mean of the time series in 
question. This set of means will then exhibit Jess scatter 
than the original series as they will be in the neighborhood 
Surrounding the point where the two regression lines 


ME and Y-f(X)] intersect the X=Y Tine. 
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4. Error Analysis 


Two sources of error, in addition to the neglected 
physics, are separation distance and the effects of friction 
in low latitudes. Investigation of the accuracy of the 
surface pressure distribution from which the wind records 
were calculated, a third possible error source, is beyond the 
scope of this thesis. 

a. Separation Distance 

The computational grid points do not exactly 
coincide with the locations of the 0S where the observations 
were made. The distance between the grid point and OS is 
defined as the separation distance. The effect of separation 
distance on the mean error, i.e., the difference between the 
mean of the observed record and that of the computed record, 
is shown plotted against separation distance in Figure 13. 

The mean error of the wind speed shows no 
relation to separation distance. However, such a relation 
is recognizable for the components of the wind vector. The 
meridional component mean error shows a slight increase with 
increasing separation distance. Although there is consider- 
able scatter, there does not seem to be any significant 
difference between Atlantic and Pacific OS in this regard. 

Mii conos not. the case for the zonal 
component however. The Atlantic OS show little mean error 
in the zonal nora At Pacific 0S, however, a rather 


pronounced dependence on separation distance is evident. 
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This dependence suggests the following: 

(1) The ten-year average of the north-south 
pressure gradient (from which the zonal wind is derived) is 
more accurately described in the North Atlantic than in the 
North Pacific Ocean. 

(2) The accuracy of this average gradient 
within 200 km of OS in the Pacific is, in part, dependent on 
the distance from the OS. 

In addition to the effect of separation 
distance on mean error, its effect on the correlation between 
computed and observed records was also investigated. The 
results, shown in Figure 14, indicate no correlation coeffi- 
cient dependence on separation distance comparable to that 
for mean error. 

b. Friction in Low Latitudes 

Low latitude wind fields play a very important 
part in large scale air-sea interaction (Berjknes, 1969; 
White and Walker, 1973). However, the wind fields in this 
region are typically difficult to model because they are 
strongly ageostrophic. As no data comparable to 0S observa- 
tions were available from low latitudes for direct verifica- 
tion, this section is included as an attempt to evaluate the 


Huscuracycot the model in this region, 
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In a recent study, Brümmer, et al. (1974) 
summarized data suggesting that the strong ageostrophic 
component of the surface wind vector in low latitudes was due 
to the effects of friction. These effects resulted in a 
reduced ratio, C, of surface wind speed to geostrophic wind 
Speed and an increase in the deflection angle, a. Table X 
lists the values of C and a as reported by Brümmer, et al., 
and those used in the model at corresponding latitudes. 

(The values of a as used in the model are for surface wind 
speeds of 2.5 m sec.) Although the values of C, as tabu- 
lated by Brümmer, et al., show no latitude dependence similar 
to those used in the model, they agree fairly well in magni- 
tude. The value of a used in the model were close to those 
tabulated by Brümmer , et al., with the exception of a ai 25. 

At this latitude the values of a and C as used 
in the model could result in an error in the orientation of 
the surface wind vector, and corresponding errors in the 
components. For example, in the case of pure zonal flow, an 
error of 30° in a would reduce the zonal component to 0.866 
Of 1tS Original value and introduce a fictitious meridional 
component equal to 0.5 of the original zonal value. Such 
Berl ctlc1ous meridional component in the computed record 
sould cad ce erroneous conclusions concerning the correlation 
of zonal and meridional wind stress anomalies. However, for 


tiem lati tudes morc Of about Sam, the parameterization of «a 


and C appear to be reasonable. 
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Br SPECTRAL ANALYSIS 

The method of analysis described in the preceding 
section allowed comparison of computed and observed values 
only as ensembles. This section describes the results of a 
comparison of frequency components constituting the observed 
and computed time series using ENEE No attempt 
will be made to discuss the atmospheric dynamics implied by 
the spectra. 

p Autocorrelation 

The autocorrelation function is the basis of 

spectral analysis. It describes the relative amplitude of 
the different frequency components present in a given time 
series (Bendat, Peirsol, op cit). The autocorrelation RE 


can be defined as: 


R. (x) s lim 


: 
; | xx (ero) y di (25) 
0 


[e 
where 


t = lag interval 
T = total time of the record 


t = time variable 


For this study T consisted of the ten years from 
1960 through 1969, corresponding to 7306 twelve hourly incre- 


ments of t. The lag interval ranged between O and one year 
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(731 twelve-hour increments). Because of the gaps in the 
observed record and the filtering performed on the time 
series (see Section II-D), the spectral analysis data are 
presented only for frequencies greater than 0.5 cycles day. 
Table XI presents, at selected frequencies, the 
values of the autocorrelation function of the computed and 
observed time series at 0SB and OSV. The frequency is deter- 
mined by the lag interval «x, which is the period of the 
frequency component in question. The discussion of the 
autocorrelation function is intended to illustrate the nature 
of the periodicity of the observed time series and the 
differences between them. (The accuracy of the model at 
specific frequencies is best described by the phase and 
coherence functions discussed in Section II-8-3.) All of the 
computed series were expected to show autocorrelation higher 
than the observed data for short lags (0 to 48 hours) due 
to the averaging nature of the geostrophic approximation 
(Munk, 1960). This, however, was not the case at either OSB 
or OSV. However, because the data were heavily smoothed at 
thes? frequencies the differences may be masked. The values 
of has C) for the observed and computed wind speed at OSB and 
OSV ere shown plotted against lag time in Figures 15 and 16. 
The annual variation in wind speed is more pronounced at OSV 
Clie menos in Spite of the fact that OSB is at a higher 
latitude. The reduced values of Rot at annual and semi- 


anneel frequencies at OSB can be attributed to a stronger 





random component in the wind speed at these frequencies. 
(0SV perhaps is influenced by the Monsoon regime as suggested 
by Malkus, 1962.) 

The nature of uy for the zonal component differs 
dramatically between OSB and OSV as shown by Figures 17 and 
18. At OSB the zonal component is essentially random for 
lags longer than 4 days. This is not the case at OSV where 
annual variation is pronounced. The flattened minimum at 
OSV (Figure 18) indicates that the ol anni periodicity 
(between 100 and 250 days) of the zonal component was reduced 
at this location by random fluctuations. 

The meridional component at both OSB and OSV, 
Figures 19 and 20, show predominately random nature for lags 
greater than 4.0 days. However, Figures 19 and 20 show tnat 
the meridional component at OSB has a slightly more periodic 
nature than at OSV. 

The computed series at both OSB and OSV show the 
Same general features as do the corresponding observed 
records. In most cases even the minor details of R(t) Tor 
the observed record are also present in the Rec) for 
computed record. This suggests that the aperiodic as well as 
the periodic nature of the wind is accurately reproduced by 
the model. 

E EECH 

For a given frequency, f, the energy density, 

G. (f) can be expressed in terms of the autocorrelation 


XX 
nction as follows: 


S 





Gf) = 2 J Ei racc (26) 


(1) » autocorrelation for lag x. This relation is 


re R 
whe D 


a Fourier transform from period space (x) to frequency space 
(f). The area under a given segment of the energy density 
curve is proportional to the variance of the time series in 
that frequency range. 

The energy density at frequencies equivalent to the 
lags (periods) listed in Table XI in the preceding section are 
listed in Table XII. A general underestimation of the energy 
by the model is indicated by this table and by the reduced 
value of the variance of the computed component records as 
Macs ted in section (111=4). Most of the "missing" variance 
in the computed records occurs at frequencies less than 


0.432 cycles day! 


(periods greater than ] week). The 
dependence of the energy density on frequency for the computed 
and observed records is shown in the upper diagrams of Figures 
2] through 26. It should be noted that the plot of the loga- 
rithm of the energy density tends to overenhance the absolute 
differences at high frequencies and to mask them at low 


frequencies. || 


The scale of the spectral energy-density plot varies 
from figure to figure as a result of attemp£ing to maximize 
mem esolutionso dates whlch varied signifücantly from plot 
o plot. 
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However, the plots indicate that the model's 
accuracy is not significantly dependent on frequency and that 
the geographical differences of the observed wind spectra 
are accurately reproduced in the computed wind spectra. 
Geographical differences are especially evident in the spectra 
of the components; this indicates that for purposes of large 
scale numerical modeling of sea-air interaction, uniform, 
zonal-averaged wind fields cannot be considered realistic. 

3x Coherence and Phase 

Coherence is a measure of the correlation between 
two time series at a given frequency; as such, it is the best 
estimator of the accuracy of the model. It can be defined 
in terms of the energy-density spectra of each of the two-time 


series of two variables, x and y, and their cross spectra. 


The equation for coherence, p is given in equation (27). 
te (f)}° 
_ y 
M S: a (27) 
S (E) GV) 


where Be) and G Wo are the energy densities of the 


y 
computed and observed time series as described in the previous 
seeRNon. the cross-spectral density function, DE is 
defined in terms of the cross correlation function between 

the variables x and y, Ee Bendat and Peirsol, EE )- 


ES Uo Uo (28) 
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Where the cross correlation function is given by: 


1 
Roy (1) = Tim J x(t)y(t*z)dt (29) 
0 
eS 
If P - ], the records are perfectly correlated at 
frequency f; for y, (f) = 0, no relation exists between x 


Xy 
and y at that frequency. 


Another estimator of the accuracy of the model 
derived from the cross-spectral density function is the 


cross-spectral phase angle, GC This function measures, 


at frequency f, the average angular difference between the 
time series of x and y, and is defined in terms of the real 
and imaginary parts of tne cross-spectral energy-aensity 


function: 


lr) tan’! E 7 0,0] (30) 


where aD: the imaginary part of Eye is called the 


quadrature spectral-density function; and bye d Ue the real 
part of Gr), is called the coincident spectral-density 


function. C and Pip can be defined as follows: 


Ee (no (31) 


AY 


tt 
NU 
ron 
C» 
sch 
Ma 
-- 
C 
— 
sch 
rr? 
== 


E AN (32) 





For this study a positive phase angle indicates 
the computed record lags the observed record in time. For 
example, a cross-spectral phase angle equal to +9° at a 
frequency of 0.1 cycles day”! indicates that the component 
of computed record at this frequency lags the observed record 
by 0.25 days or 6 hours. 

The values of the coherence and phase (cross- 
spectral phase angle) at OSB and OSV for the frequencies of 
interest are tabulated in Tables XIII and XIV. In Figures 
21 through 26 the coherence and phase are plotted in the 
lower two diagrams. At both OSB and OSV there is a gradual 
increase in coherence with decreasing frequency. AS Suggested 
in Section III-A-3, the coherence between computed and 
observed wina speed is lower than is the coherence between 
the observed and computed components. At OSB there is a 
constant increase in Dt with increasing frequency. For 
the wind speed and the components of the wind vector this 
phase difference is approximately 45° at f=0.5 cycles day !, 
indicating the computed record lags the observed by about 6 
hours. 

For OSV there is no consistent phase difference 
between the observed and computed records similar to that at 
OSB; however, the coherence is similar to that for OSB, and 


approaches unity at periods of six months and a year. 


29 





The data presented in this section suggest that 
this model reproduces the major features of the observed 
energy density spectra of wind speed and orthogonal components 
for frequencies between 0.5 and 0.0027 cycles day"! (periods 
of 12 hours to one year). It also suggests that spectra for 
the components of the wind vector T be deduced from the 


corresponding spectra of the wind speed. 
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IV. SUMMARY 


A. ACCURACY OF THE MODEL 

The results of this study suggest the following: 

ES The relation between the computed and observed 
records of wind speed and components of the wind vector is 
essentially linear. 

2. The model exaggerates both high and low values of 
wind speed and components and this exaggeration is most 
pronounced at low latitudes. 

4 The accuracy of the model as shown by correlation 
with observations is slightly higher in the Atlantic compared 
to the Pacific Ocean; and it is. higher in mid-latitudes than 
in either high or low latitudes. 

4. The ten-year (1960-1969) trend of wind conditions 
between 29°N and 56.5°N is similar for both the computed 
records and the observed records; and observations indicate 
a slight reduction (-0.5 m sec!) in wind speed over the 
northern hemisphere between 35N and 50N. 

3 The model explains between 30% and 70% of the 
variance of the parameters measured at 12-hour intervals. 

6 There is a fictitious southerly component incorpo- 
rated in the computed records at all grid points, possibly 
the Esult of a fictitious east-west gradient in the pressure 


analyses used. 
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T: The deflection angle between the surface wind and 
the geostrophic at latitudes between 0° and 5°N as used in 
the model is substantially less than the values thus far 
observed. 

Se The accuracy of mean wind conditions derived from 
the model increases with the length of the averaging interval, 
as indicated by increasing coherence with longer period 


fluctuations. 


B. POSSIBLE APPLICATIONS 

The last conclusion points to a possible application of 
the model in addition to its current operational use. Monthly 
means of wind conditions derived from 12-hourly surface 
pressure analysis using the model should exhibit even greater 
accuracy than do the 12-hourly values. The mean wind para- 
meters derived in this way would be true averages and not 
resultant wind parameters. A historical record of this 
nature should permit a more realistic quantitative approach 


to sea-air interaction studies than has been the case to date. 
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TABLE I 


Monthly distribution of surface pressure analyses. The 
number of analyses used is in the numerator, and the 


denominator equals the number of twelve hour intervals in 
the month. i 
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TABLE IX 


Scatter diagram symbols (for Figures 11 and 12) 
and the corresponding number of coincident 
computed and observed wind speed pairs. 


Diagram 


Symbo] 
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JD Ob op ma Evo os OU GA bh 


Equivalent 
Values 
al J 20 
2 K 21 
3 L 22 
4 M 23 
5 N 24 
6 O 26 
7 P 27 
8 Q 29 
9 R 33 
10 S 41 
11 T 57 
12 U 89 
13 V 159 
14 W 2909 
15 A 59g 
16 Y 1049 * 
17 Z 2073 
18 * Al 
19 
TABLE X 


Ratio of surface wind speed (Vs) to geostrophic 
) and deflection angle (in degrees) 
des as tabulated by Brümmer, et al. 
in the model. 


wind speed (V 
at low latitu 
and as formulated 
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TABLE XIII 


- 


Coherence between observed and computed time series 
for Ocean Station Bravo and Ocean Station Victor 
at selected frequencies. 
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TABLE XIV 


Phase difference (in degrees) between observed and 
computed time series for Ocean Station Bravo end 
Ocean Station Victor at selected frequencies. 
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Figure 1. Latitude dependence of the ratio of the 
surface wind speed (Vs) to the geostrophic wind 
speed (Vg) as found by Roll (broken line), Carstensen 
B line) and às formulated in the model (dashed 
line 
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Figure Z. Latitude dependence of the d 
for various wind speeds (M in m sec-! 
in the model. 
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Figure 3. Latitude dependence of the deflection angle 
cuna py Roll (broken line), Carstensen (solid 
line). anamas formulated in the model for Wind speeds 
mono. m sec-| and 15 m cec l (dashed lines). 
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Etre 4. Result of five point binomial filter applied 
to a typical portion of a zonal component time series. 
Original series (solid line) and filtered series 
Usted line) are shown in the upper graph. The 
residual series is shown in the lower graph. 
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Pigure 8. Total change (integrated linear trend) 
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December 1969 in the observed (solid line) and 
computed (dashed line) wind speed and zonal and 
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eneede (in msec !) at Ocean Station Wevember. 
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Figure 21. Log energy density (in m^ sec 
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wind speeds at Ocean Station Victor. 
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